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ABSTRACT 

Radiative accelerations on Ne are calculated for the atmospheres of main sequence 
stars with 11000 < Toft < 15 000 K. This range corresponds to that of the mercury- 
manganese (HgMn) stars. The calculations take into account neon fine structure as 
well as shadowing of neon lines using the entire Kurucz line list, bound-bound, bound- 
free and free-free opacity of H, He and C as well as some NLTE effects. NLTE effects 
are found to modify the radiative acceleration by a factor of order 10^ in the outer 
atmosphere and are crucial for dm < gcm~^. The dependence of the radiative 

accelerations on the Ne abundance, effective temperature and gravity is studied. Ra- 
diative accelerations are found to be well below the gravity in the entire range of Tcs 
and it is predicted that in stable atmospheres devoid of disturbing motions, Ne should 
sink and be observed as underabundant. This agrees with recent observations of low 
Ne abundances in HgMn stars. 
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can be followed in time on l y in the en ve lope, as shown in 
the case of neon by ^eaton ( 1997 , 199S: ); Landstreet , Dole z 



Vauclair| ( |l998| ) and [LeBlanc, Michaud fc Richeij ( pOOO^ 



upper main sequence CP stars (Michaud 1970). The up 



ward radiative acceleration and downward gravity acting on 
different chemical species compete and drive their micro- 
scopic diffusion, resulting in the stratification of these ele- 
ments in the stellar atmospheres or envelopes of these stars. 
If these atmospheres are sufficiently stable, chemical inho- 
mogeneities can persevere, and are not wiped out by various 
mixing processes such as convective zones, stellar winds or 
rotationally induced mixing. In the particular case of HgMn 
star atmospheres, temperatures are high enough for hydro- 
gen to be largely ionized, thus removing its convection zone. 
Also, if the He ionization convection zone diminishes due to 
He settling in la t e B stars rotating slower than c. 75 kms^^ 



These calculations suggest that radiative accelerations on 
Ne are much lower than gravity below the atmosphere and 
Ne should sink and be observed as underabundant unless a 
weak stellar wind transports upwards the Ne accumulated 
in deep layers. 

Nevertheless, observations are restricted to the atmo- 
spheres. This important stellar region bridges the effects in 
the envelopes with those outside the stars, such as stellar 
winds, and one clearly needs also to know the radiative ac- 
celerations here in order to compare the theory with observa- 
tions. The case of Ne is worth studying for two reasons: it is 
an important element having very high standard abundance 



(Michaud 1982; Charbonneau fc Michaud 1988), diffusion 



([Anders fc Grevesse| |l989 ; Grevesse, Noels & Sauval 



199e 



Dworetsky fc Budaj| pOOC ) of about A(Ne/H) = 1.23 x 10" 



can operate very effectively in the atmosphere. This is the 
case for HgMn stars as they are also known as slow rotators. 
Such diffusion processes are generally time dependent prob- 
lems, but the characteristic time-scales in the atmosphere 
are much shorter than in the envelope (< 10^ yr) due to 
the considerable lower density, smaller collision rates and 
higher diffusion velocities (c. 1 cms~^). This makes the time 
dependent calculations extremely difficult. So far, diffusion 
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and because one can expect similar effects as in the case of 
He. Nel (along with Hel) has a very high ionization poten- 
tial and first excited level well above the ground state. Ne I 
is thus the dominant ionization state in the atmospheres of 
these stars and all Nel resonance lines are in the Lyman 
continuum. This part of the spectrum is very sensitive to 
non-LTE (NLTE) effects and the temperature structure of 
the atmosphere. It is known that there is a big difference in 
the temperature behaviour between LTE and NLTE model 
atmospheres at small optical depths and that temperature 
can rise outwards in the NLTE models. Nel is also known 
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to have strong NLTE effects on lines originating from ex- 
cited levels in the visible part of the spectrum even i n late 



B-type main sequence stars (Auer fc Mihalas 



199£ 



197S 



Sigut 



Dworetsky & Budaj 2000p . Thus, the main compli- 



cation is clear in doing the task. One cannot rely upon the 
diffusion approximation when calculating energy ffux as this 
is a non-local quantity in the atmosphere, in contrast to en- 
velope calcula tions. One needs to solve the radiative trans - 
fer as, e.g., in Hui-Bon-Hoa, L e Blanc & Hauschildt _[ 2000 



1996 



Hui-Bon-Hoa, Alecian fc Artru 1 
( ^002| ), including NLTE effects if possible 



Hui-Bon-Hoa et al 



Unfortunately, observations of Ne in HgMn stars are 
very scarce and apart from rec ent LTE Ne abundances in 
K Cr ic, HR7245 and HD296 47 by |Adelman fc Pintado| (|200C| ) 



and Adelman et aL ( 2001) the only study of Ne was under 
taken recently by [Dworetsky fc Budaj| ( |2000[ ). Their NLTE 
analysis of 20 HgMn and 5 normal stars revealed that Ne is 
clearly underabundant in HgMn stars and the deficit reaches 
its peak values in the middle of the T^s domain of HgMn 
stars. 

In this paper we deal with the radiative accelerations of 
neon in the atmosphere. This might enable us to model the 
vertical element stratification in the future, e.g. considering 
the atmosphere as being in a sequence of stationary states 
which are to be tailored to the much more slowly varying 
envelope as proposed by Budaj et al, (1993). Nevertheless, 



some conclusions can be drawn already from the radiative 
accelerations alone. In the following, cgs units and LTE ap- 
proximations (in a NLTE atmosphere model) are used if not 
specified otherwise. 



2 RADIATIVE ACCELERATION 
2.1 Theory 

The physics and calculations of ra diative acceler ations in the 

(IioqtI), 



envelopes are clearly explai ned in 



provements are described in Gonzales et al 



Seaton 



theoretical aspects are discussed in Bapar & Aret (1995) 



(|1995 ) and some 



Our approach to the problem of calculating accelerations in 
the atmospheres is described here. Radiative acceleration on 
the element is due to the momentum transferred to the ele- 
ment ions from photons absorbed via bound-bound, bound- 
free and free-free transitions. Bound-bound transitions are 
usually the most important. An expression for the radiative 
acceleration on the ion i through its spectral lines can be de- 
rived from the radiative transfer equation. For a particular 
spectral line in the case of complete redistribution we have: 



dl,y = [nu,iAui — {nijBiu — nu,iBui)Iv]hvipiu{v)ds 



(1) 



where the three terms describe spontaneous isotropic emis- 
sion, absorption, and stimulated emission, respectively. Note 
that stimulated emission is in the direction of the photon 
absorbed. Here, and in the following, Ii, is monochromatic 
intensity, Fi, is monochromatic energy flux, is the energy 
of the transition from the lower level I to the upper level u, 
(fiiu{v,T) is the normalized Voigt profile, r is the Rosseland 
optical depth, s is the distance along the beam, ni^i is the 
NLTE population of the l-ih state of i-th ion, rii is the total 
NLTE i-th ion population, bi^i,bu,i are b-factors defined as a 
ratio of NLTE to LTE level population {bi,i = ni^i/n*^^), m is 



the mass of the i-th ion, c is the speed of light, k is the Boltz- 
mann constant, and T is the temperature. Aui, Bui, Biu are 
the Einstein coefficients where Bi^ is related to the oscillator 
strength fiu and Bui by: 



Biu = vre fiu{mcchu) 



9l,iBlu — gu,iBui 



(2) 



where e, rric are the electron charge and mass respec- 
tively and gi^i,gu,i are statistical weights of the lower and 
upper levels. Introducing the vertical geometrical depth 
dr = cos 6ds, multiplying Eq. ^ by cos ddui and integrating 
through all solid angles duj we get for the radiation pressure 
gradient dpd/dr: 



dpi 
dr 



fT'U,i Ql^i 



ni^iBiuhvF„ifiu{v) 



(3) 



If Pill is the momentum removed from the radiation field 
through one spectral line in unit volume and aiu, a,i are the 
corresponding accelerations acquired by the ion i then we 
get: 



E 

lu 



1 v-^ dpiu 
rarii ^ — ' dt 

lu 



ai{T) = — / '^^dv 

mrii ^ — ' /„ dr 



(4) 



(5) 



The sum runs through alH ^ ti transitions in the ion. Ap- 
plying the Boltzmann formula to LTE populations, in this 
general NLTE case, we finally have: 



^ mc 



lu 



6 I ^lu~ 



Fv(T)ipiu{v,T)dv 



(6) 



Note that, the particle fiux of the element associated with 
such definition of radiative acceleration would be: Q 

(7) 



J = ^ n-iVi = ^ riiDi^Ui 



where Vi,Di are the i-th ion diffusion velocity and diffusion 
coefficient, respectively. 

The above expression for radiative acceleration includes 
the correction for stimulated emission w hich is o ften omit- 
ted by other authors, as pointed out by Budaj (1994) and 
Seaton (1997), b ut does not include the r e distri bution effect 
pointed out by Montmerle fc Michaud ( 197(:) a r id ge ner- 
alized to include gravity by lAlccian fc Vauclaiil (|l983|). I t 
has been substantially revisited by Gonzales ot al, (1995). 
The main idea of the redistribution effect is that, after the 
photon absorption, ion i is in the excited state u and has 
non-negligible probability (1 — r„,i) of being ionized (mainly 
by collisions with fast electrons) before losing its momen- 
tum (mainly by collisions with protons). The probability of 
its remaining in the state i the probability of re- 

combination or further ionization is negligible. Momentum 
absorbed during the transition in the state i, piu, should 
then be redistributed to the two mass reservoirs mni+i and 

^ O ne can d efine the radiative acceleration in some other way as 
e.g. [jeaton| (1997), corresponding to different expressions for the 
particle flux. 
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mrii of 4 + 1 and i ions in the proportions (1 — ru,i) and r„,i, 
respectively. Consequently, one gets: Q 



mrii 



lu 



dt 



E 

lu 



(8) 



mjii+i ^-^ at ^-^ rii+i 

l\i lu 



(9) 



The redistribution function, r„,i, depends on the ionization 
rate from the particular level, l3u,i, and the collision rate of 
the ion, /3i: 

A 1 



l3u,i + Pi 



+ 1 



(10) 



It is often useful to picture the radiative accelerations 
on many ions being applied via some effective radiative ac- 
celeration on the element as a whole. Consequently, to com- 
pare the results with other authors and to display them, 
we also calculated the weighted mean radiative acceleration 
(although the results are stored for each ion and abundance 
sep arately). We adopted the following expression suggested 
by Gonzales et al. (1995): 



a(r) 



(11) 



Di is obtained from -Di,Hi, -Di.hii, which are diffusion coeffi- 
cients which describe the diffusion of the i-th ion in neutral 
and ionized hydrogen, respectively: 

1 _ 1 1 



A, 



A, 



(12) 



For the diffusi on of the i-th ion (exclu ding neutral) in pro- 
tons we used (AUer & Chapman 1960): 



D 



1.947 X lO'T 



9rp5/2 



i,mi 



(13) 



nmiM2ZfAi{2) 

where nnii is the proton number density, Zi is the ionization 
degree (0 = neutral, ...) and: 

A 



yli(2) =ln(l + x^,) 



M, = 



1 + A 

where A is atomic weight in atomic mass units and: 
4dDfcT 



Xu 



and du = 



kT 



where do is Debye length and is electron numb er density. 



For the diffusion of n eutrals in protons we used (Michaud 



Martel fe Ratel 



197J 



A.Hii = 3.3 X 10^ 



riHii 



1 + A 
aA 



(14) 



where a is the polarizability of neutral species {a = 
0.395,0.667 x 10~^'*_cmf for Nel and Hi respectively, from 



Teachout & Pack (1971)). The same formula was used for 
the diffusion of ions in neutral hydrogen with a a(HI) 

^ Note that these formulae are valid for our radiative accelera- 
tions, i.e. as defined by Eq. ^. Other formulae for the redistribu- 
tion effect without factor would also hold but if combined 
with a sli g-htlv differ ent definition of radiative acceleration, e.g. 
that of Beaton (1997). 



and TiHii nm, where nm is the neutral hydrogen number 
densit y. For the diffusion of n eutrals in neutral hydrogen we 
used (Landstreet et al, 1998): 



5.44 X lO^VT 1+A 
nm{Sm + S,)''\l A 



(15) 



where 5ui,Si are If I and diffusing atom diameters, respec- 
ti vely. The following values calcula t ed us ing the method 
of Clementi, Raimondi & Reinhardt ( 1963[ )P] were adopted: 



1.06 10^*^ cm, 5noi = 0.76 10^^ 



2.2 Line selection and numerical calculations 

As one goes from the envelope higher into the atmosphere, 
the electron density and associated Stark broadening drop 
sharply and the spectral lines become much narrower. One 
needs to consider fine structure to integrate properly the 
contributions from bound-bound transitions. Atomic data 
for t he Nel-iv lin es were extracted from the Kurucz line 
list ( |Kurucz|p^90| , CDROM 23). However, one usually does 
not need to consider all the lines of the element and it is of 
high practical importance to select just the most important 
ones. Lines which are important for a particular atmosphere 
model were selected in the following way. 

For practical purposes it is convenient to ascribe three 
characteristic Rosseland optical depths ti,T2,T3 or charac- 
teristic temperatures Ti(ri) < T'2(t"2) < 73(t"3) to each ion 
i of the element of interest. The range Ti to T3 spans the 
region where the ion i is sufficiently populated and its con- 
tribution to the total radiative acceleration of the element 
is thus not negligible, while T2 is the approximate tempera- 
ture at which the ion dominates the element and the i-th ion 
opacity is maximal. For each line lu of the ion i two parame- 
ters RAiju.i and RAi^i^^s were calculated which correspond 
to two characteristic temperatures Ti , T3 : 



RA,juA,3 = 8.853 X 10' 



' [1- 



(16) 



where Eij is the excitation potential of the lower level and: 



F4ti,3) = 7rB,(Teff) if Ti,3 < T.g, or 



4 ^ dBAT)dT 



if Ti.a > Tef 



(17) 
(18) 



3 Xc,,. dT dr 

where Bi, is the Planck function, T^b is the effective tem- 
perature, and Yc.i/ the approximate continuous opacity after 
Borsenberger, Michaud & Praderie (1979). 



The parameter 1,3 is essentially a rough estimate 

of the expected radiative acceleration through weak lines, as 
Eq.^ takes only opacity in the continuum into account and 
partition functions are set to 1.^ Then we calculate a sum 

^ http:/ /www. webelements.com/webelements. html 
^ The first simplification does not have any serious effect on this 
selection of the most important lines of a particular ion as it 
means we overestimate the expected accelerations through strong 
lines, so we cannot omit them. The second simplification also 
does not have any effect on line selection of an ion as RAi 1 3 
or radiative accelerations for all the lines of a particular ion are 
scaled by the same value of ion population which is proportional 
to its partition function. Again, it also means we cannot omit 
important lines. 
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through all the lines of the i-th 



16000 



ion and skip the lines which do not contribute significantly 
to the sum. We ended up with about 10^ Ne lines for a 
particular model atmosphere. 

To evaluate numerically the integral in Eq. ^ it is neces- 
sary to estimate the integration step and integration bound- 
aries. Thermal Doppler broadening at the coolest layer and 
classical radiative damping put constraints on the width of 
any spectral line and we use minimum frequency spacings 
equal to 



& 12000 



^ (Ai.D(ri) +rR/47r)/2 



(19) 



8000 



where: Ai^d ~ vc~^ ^j2kTi/m and Fr 
Integration boundaries are associated with our line width, 
which is generally the broadest at the place of highest elec- 
tron density due to Stark broadening and we use: 



T.,^12000 K, log E=4 



Paschen cont. 



Lyman cont. 



NLTE 



10'° 



LTE^ 



W lo- 

dm [cgs] 



10"" 



10"' 



Ai^w = fci[r(r3)/47r + Ai^D(r3)] 
where 

r(r3) = rR + rs(r3)-i-rw(r3) 



(20) 



(21) 



is the sum of the radiative, Stark and Van der Waals damp- 
ing parameters. Here, k\ is an adjustable constant for which 
numerical tests revealed that fci « 2 — 3 is high enough. This 
is, however, just the case of a weak line. If the opacity in the 
center of our line is much greater than the opacity in the 
continuum one must usually integrate much further until the 
line opacity drops significantly below continuum. The reason 
why is apparent if one expresses the flux in Eq.^ following 
the difi'usion approximation and neglects continuum opac- 
ity; it is inversely proportional to the Voigt function. Conse- 



Figure 1. Temperature behaviour in a representative Kurucz's 
line blanketed LTE and Hubeny's NLTE model atmosphere. Ver- 
tical bars indicate the unit monochromatic optical depth in the 
Lyman and Paschen continua close to the important Nel A736 
and Ne I A6402 transitions. As a depth coordinate we use 'dm' (g 
cm^'^) which is the mass of a column above a unit area at that 
depth in the atmosphere. 



of Ne have been performed using A6402 (Dworetsky & Bu- 



quency. Only when the integration is carried far enough into 
the line wing, where continuum opacity again dominates, 
does flux approach a constant function of frequency and the 
integrand sharply fall. Finally, the integration boundaries 
should be the maximum value of the mentioned effects: 



quently ~he integrand in Eq.^ is a constant function of fre- daj 200Cl| ). The Lyman continuum is very opaque and models 



which go very high in the atmosphere are needed to solve 
properly the radiative transfer at these wavelengths. 



Fu{t) in Eg. o was calculated usin 



g the SYNSPEC36 and 



Ai^^ 



±MAX[Ai.w(r3);Ai.,(r2) 



(22) 



SYNSPEC42 ( Hubeny, Lanz fc Jeffery 1994 ) codes to solve 
the radiative transfer modified to get a nonstandard out- 
put of flux at each frequency and depth (I. Hubeny, priv. 
comm.), and to take into account the elements with atomic 



where Avs{t2) is defined for strong lines only and is for 
which the line opacity drops significantly below the contin- 
uum opacity at the depth where the involved ion population 
reaches its peak values. Typically, anything from 40 to 10^ 
frequency points were necessary to calculate the radiative 
acceleration through a Ne line. 

The models adopted here wer e NLTE mod e l atmo- 



spheres calcu l ated with tlusty195 ( [Hubeny 
& Lanz 



1992 



1988 



Hubeny 



1995) for standard solar composition and zero 



microturbulence. We treated H I, He I, He ll, C I, C ll, C ill as 
explicit ions, which means that their level populations were 
solved in NLTE under the assumption of statistical equi- 
librium and their bound-bound, bound-free, and free-free 
opacities were taken into account. Carbon was included as 
it is an important opa city source at cooler effective temper- 
atures (|Hubeny||l98l|). In Fig. |l| we plot the NLTE model 
against the [Kurucz (199S) LTE line blanketed model for 
the same effective temperature and gravity. The two verti- 
cal bars point to the area where the continua around Nel 
A736 and Nei A6402 lines are formed. The first is a strong 
resonance line in the Lyman continuum and the second is in 
the Paschen continuum in the visible. Abundance analyses 



number > 30 (Krticka 1998). This enabled us to take into 
account the bound-free and free-free opacity from explicit 
levels specified in the atmosph ere mod el and the line blend- 
ing using the full line list of Kurucz (1990). All the Ku- 
rucz lines with line-to-continuum opacity ratio greater than 
10'** at r = MAX[1 , T2] were considered. All the resonance 
lines of Ne and other elements were calculated in approxi- 
mate NLTE (allowing the source function to deviate from 
B y followi n g the 'second-order escape probability method' 
of |Rybicki| ( |l984| ), see Eq. ||). 

Radiative accelerations were then calculated following 
Eq. ^ assuming bu,i ~ bi^i — 1 and LTE level popula- 
t ions. Partition functi ons were taken from the uclsyn code 
( ^mith fc Dworetsky 



1988 



Smith 1992). Radiative accel- 



erations were calculated for Nel-iv and four abundances: 
50 X Aq, Aq, IQ-'^ X Aq, 10'^ Aq where Aq = 1.23 x 
lO"** and for the following model atmosphere parameters: 
Teft = 11 000, 12 000, 13 000, 14 000, 15 000 K and log g = 4, 
plus models with = 12 000 K and log q = 3.5, 4.5. These 
are av ailable in digital form in Budaj, Dworetsky & Smalley 
(|oo|) including a fortran77 code containing the partition 
function routines. 
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Figure 2. Nc ionization fractions in LTE in a representative 
NLTE model atmosphere. 



3 RESULTS AND DISCUSSION 

3.1 Radiative accelerations on different Ne ions 

The LTE populations of various Ne ions throughout the at- 
mosphere are depicted in Fig. ^ for a representative atmo- 
sphere model with Teff = 12 000 K, logy = 4. Nel is defi- 
nitely dominant in the atmosphere. Ne iv becomes populated 
at the base of our models while Ne ll may dominate at very 
low densities and electron concentrations for dm < 10~* 
gcm~^. But above dm < 10~^gcm~'^ strong departures 
of N e ionization and excitation equilibrium from LTE oc- 
cur (Sigut 1999) and our accelerations on different Ne ions 
should be considered as approximate above this level (see 
Fig.^ for an illustration of expected deviations from fully 
consistent NLTE case). 

Accelerations in the representative model are plotted in 
Fig. ^. One can see that none of them exceeds the gravi- 
tational acceleration except at the very deepest layers. The 
acceleration on Ne I is strongest in the atmosphere where this 
ion is also most populated. The contribution of the neutral 
species like Nei to the total effective radiative acceleration 
of the element is both favoured by weighting with their dif- 
fusion coefflcient (see Eq. |ll|) which is approximately two 
orders of magnitude larger than that for the charged species 
and reduced by the redistribution effect. Nevertheless, ac- 
celeration of Ne I will be the most important contribution to 
the total acceleration not only in the atmosphere but also at 
rather great depth up to dm ~ lOgcm"^. The low radiative 
accelerations found on Nel are a consequence of its atomic 
structure because all resonance transitions, which are usu- 
ally the most important, are in the Lyman continuum where 
Nel atoms see little photon flux, while the rest of the lines 
originate from highly excited and weakly populated states. 
Nor can the redistribution effect discussed earlier increase 
the Ne I acceleration (see Eq. ^ and Fig. just the oppo- 
site occurs. Radiative acceleration gained in the Nel state is 
redistributed between Nel and Nell, which reduces the Nel 
acceleration and adds it to the Nell state where, however, 
the acceleration is not so effective due to the much lower 
diffusion coefflcient of ionized species; see Eqs. |^, pj] . 

We would like to stress that radiative acceleration (Eq. 
p|) on a particular ion is very sensitive to the partition func- 




> 



10"^ 10"' 10"' 
dm [cgs | 

Figure 3. Radiative accelerations on Nel-iv ions in a represen- 
tative NLTE model atmosphere for the standard abundance of 
Ne. The downward gravitational acceleration of 10** cms"'^ is at 
least 2-3 orders of magnitude greater than the upward radiative 
acceleration throughout the photosphere. 



tions. These are generally precise in the region where the ion 
dominates the element but have large uncertainties outside 
the region. Displaying effective accelerations overcomes the 
problem as partition functions drop out in Eq. ^ outside the 
region. Thus, if the reader wishes to use our radiative accel- 
erations stored separately for each ion, the same partition 
functions should be used also. 



3.2 Dependence on effective temperature and 
surface gravity 

Now we can proceed further and explore how the situation 
changes with the effective temperature of the star. It is clear 
from Fig. |^ that the total acceleration rapidly increases with 
the effective temperature. However, within our range of in- 
terest it never exceeds the gravity in the line forming re- 
gion. Nevertheless, this suggests that Ne could accumulate 
at deeper layers in cooler stars than in hotter ones, the lat- 
ter thus being more vulnerable to departures from the ideal 
stable atmosphere. Various mixing processes or radiatively 
driven stellar winds w hose intensity also increases with th e 

[Krticka fc Kubat| pOOl]) 



1995 



effective temperature (Babel 
might more easily enrich Ne in the atmospheres of much 
hotter stars. This resembles the case of He which has atomic 
structure similar to Ne. The abundance of He progressively 
increases with effective temperature from He-weak towards 
He-rich stars. This suggests that it might be interesting to 
search for Ne-rich analogues of He-rich stars among early 
B stars. On the other hand, for Toff < 11 000 K, hydrogen 
becomes partially ionized, which triggers ambipolar diffu- 
sion and hydrogen superficial convective zones. This may 
also transport some Ne to the atmospher e from below. Our 
observations ( pworetsky fc Budaj 2000) indicate that Ne 
underabundances are most pronounced in the middle of the 
HgMn domain and tend to be less extreme towards the cool 
and hot ends of the HgMn temperature region. This pattern 
could therefore be qualitatively expected. Observing Ne be- 
low Tcff = 11 000 K, however, becomes very difficult. 

The radiative accelerations also depend on the surface 
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10"* 10"'' 10"** 10"- 1 10- 

dm [cgs | 



Figure 4. Radiative accelerations on Ne throughout the atmo- 
sphere as a function of stellar effective temperature and depth. 
Calculated for standard abundance of Ne. 



1 1 1 1 


10"5 // ■ 


T^fP 12000 K, log g=4 




A(Ne) = (50, 1, 10"^ 10"^) x A^,j(Ne) 




\ / 





10"** lO"*" 10"** 10"^ 1 10^ 

dm [cgsj 



Figure 6. Radiative accelerations on No throughout the represen- 
tative NLTE model atmosphere as a function of the homogeneous 
Ne abundance. 



=12000 K 



log g=3.5, 4.0, 4.5 




10"° 10"" 10"' 10"- 1 10- 

dm [cgs] 

Figure 5. Radiative accelerations on Ne throughout the atmo- 
sphere as a function of stellar surface gravity and depth. Calcu- 
lated for standard abundance of Ne. 



gravity (Fig. |^. Lowering the surface gravity of the model 
shifts the radiative acceleration behaviour to deeper layers 
and reduces the radiative acceleration at the bottom of the 
atmosphere. However, the gravity to radiative acceleration 
ratio does not change very much and in this context effec- 
tive temperature is a more important atmospheric param- 
eter than log jr. Apparently, within the main sequence (ap- 
proximately 3.5 < logg < 4.5), radiative acceleration of Ne 
never overcomes gravity. 

3.3 Effects of homogeneous Ne abundance 

Radiative acceleration on the element also depends on its 
abundance. This is a simple consequence of the fact that 
the energy flux in the integral in Eq. |^ depends on the el- 
ement abundance. This is, however, only important in the 
case of strong lines where the line opacity is comparable to or 



larger t han opacity in the continuum. One can then expect 
that radiative acceleration will asymptotically increase as el- 
ement abundance decreases, until it reaches some maximum 
value corresponding to the sufficiently low abundance when 



all the element lines disappear and radiative acceleration 
will no longer be sensitive to t he abundance (for more detail 

The situation is illustrated 



see Alecian 



LeBlanc 



200C) 



in Fig. ^ Surprisingly, the accelerations are very little de- 
pendent on abundance in the photospheres above dm ~ 0.1 
gcm~^ and do not exceed gravity for any abundance except 
at the extreme base of our models. Here a Ne deficit rang- 
ing from -0.5 dex for the hottest to -2 dex for the coolest 
model could be supported by radiation (i.e. radiative accel- 
eration roughly equals gravity for this abundance and depth, 
dm ~ 20gcm~^). This weak abundance sensitivity results 
from the fact that resonance Nel lines are not very strong 
because they are in the Lyman continuum where the contin- 
uous opacity is very large, effectively reducing line to opacity 
ratio. We have found that contributions to the radiative ac- 
celeration of lines originating from excited Ne I levels is also 
important, but these weak lines do not introduce strong de- 
pendence on the abundance because their line-to-continuum 
opacity ratio is low due to the high excitation energies and 
low populations of the levels from which they originate. One 
can see from the figure that there is no abundance for which 
the radiative acceleration could balance the gravity in the 
line forming region above dm, ~ 0.1 gcm~^. The implication 
is that in this region, assuming a stable atmosphere devoid 
of any motions, Ne should sink and should be almost com- 
pletely depleted and only its strong concentration gradient 
could balance the downward flux of Ne atoms. The Ne atom 
is at least an order of magnitude heavier than the mean 
'molecular' mass of the gas, so the characteristic scale of 
the Ne abundance gradient would be an order of magnitude 
less than the pressure scale height, i.e. log dm ~ 0.1. Con- 
s equently, the fact that Ne is detected in most HgMn stars 
(Dworetsky & Budaj 2000) suggests the presence of some 



kind of instability competing with diffusion and suppressing 
the efficiency of diffusion processes. An example of such an 
instability is the weak stellar wind suggested by Landstreet 
et alj (|l998|). 



One may also notice a sudden drop in the radiative ac- 
celerations at approximately 10~^ < dm < 10~^ gcm~^. 
This is not a numerical artifact but a very interesting ef- 
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Figure 7. An illustration of the redistribution effect on the ra- 
diative accelerations on Ne throughout the representative NLTE 
model atmosphere as a function of depth calculated for standard 
homogeneous abundance of Ne. (a) — no redistribution - solid 
line, (b) — no r edistribution and no d iffusion coefficient weighting 
— dashes, (c) — Gonzales et al. ( 1995| ) type redistribution - dots, 
(d)- 



type rcdiotribution — daoh doto. The i 



' pointo to the 



depth ..^herp T = ■■^innnA' 



feet when radiative acceleration may acquire negative val- 
ues; this will be discussed in more detail below. 



3.4 Redistribution effect 

To take the redistribution effect into account properly one 
would need to calculate probabilities of all possible radia- 
tive and coUisional transitions from the upper level to con- 
tin uum including multipl e cascade transitions as suggested 



by Gonzales et al, (L995). That is a very difficult task and 
the above authors developed a method in which they set 
ru,i ~ for n > 3 and ru,i ~ 1 for lower states. Thus 
their r„,i is a simple step function of level energy and does 
not depend on e.g. temperature. To cope with the prob- 
lem and estimate the influence of the redistribution effect 
we propose a slightly different and more gen eral approach. 
The idea is simple. As mentioned in Section 



2.1 



collisions 

with electrons may change the ionization state before the 
momentum absorbed by the atom is lost in collisions with 
protons. Consequently, the redistribution function from a 
particular excited level depends mainly on the electron to 
proton coUisional rates ratio. Higher atomic levels with en- 
ergies corresponding to the electron kinetic energy ~ kT 
below the ionization threshold are most strongly coupled 
to the continuum via coUisional excitation and ionization. 
Probabilities of such processes are strongly te mperature de- 



1962 



pendent via the Boltz mann factor e"^^/*^^ (Seaton 
Van Regemorter 1962|) and still deeper atomic levels are af 
fected as one goes deeper into the star. Electron to proton 
number density ratio almost does not change with depth in 
such hot stars and the Boltzmann factor thus embraces the 
essence of the electron to proton collision rates ratio and 
one may write Pu.i/Pi = Ce"^^ ^ ^'^'^ . Consequently, based 
on Eq. holwe suggest the following expression for the r^.i- 



Ce-u + 1 



with U = {li- E, 



1 



BkT 



where li is the i-th ion ionization potential and A, B, C are 
adjustable constants. It is beyond the scope of the present 
paper to calibrate precise values for the above mentioned 
constants. B is most important as it parameterises the tem- 
perature dependence and is of the order of 1. The 1/A term 



was in troduced for those who prefer to use a Gonzales et al 
( 1995| ) type of redistribution but do not like its step function 
shape and A is of the order of kT, but we use A — > cx3. C 
parameterises the ratio of the total electron to proton col- 
lision rates. If Eu,i h the level will register all electron 
impacts and I3u,i/I3i = C. Because electrons are -y/ rnp/me 
faster than protons, their collisions will be more frequent by 
approximately the sa me factor, and C is thu s of the order 
of 40. The method of Gonzales et al, (1995) is equivalent 
to the approximation of Eq.|23| by a temperature indepen- 
dent rectangle and could be approached as a special case 
when B — > oo and A — > or A — > oo depending on whether 
Eu,i > li — Ecrit or Eu,i < h — Ecru, respectively, where 
Ecrit is some threshold energy. Note that ^4 — > means 
ru.i 1 i.e. no redistribution at all and A —> oo means 
ru,i ^ 1/(C' -I- 1) — > a small value, i.e. total redistribution is 



approach e d bu t never fully realised. Assuming that Glonza- 
les et al. (1995) took the lowest n = 3 level of C ll as a cut- 
off in their calculation for T = 31000K which corresponds to 
Ecrit = 80000 cm~^ below continuum, we can simulate their 
method (using our Eq. ^ ^ and forcing ru,i = 1 or ru,i ~ 0) 
and calibrate our method so that it would give ru,i ~ 0.5 
at the same temperature. Assuming C — 40, A ~* oo we get 
from Eq. |2^ that B = 1.0 - a reasonable value indeed. Re- 
sults of both methods are illustrated in Fig. |^ and one can 
see that while bo th methods may give similar results for 
T = 31000K the IConzales et 



al 



(1995) method strongly 



overestimates the redistribution effect for lower tempera- 
tures. A simple test of the redistribution effect treatment 
is recommended by calculating radiative accelerations with- 
out weighting them by diffusion coefficients (setting Di — 1 
for example in Eq. ^l|). This is because such accelerations 
are an invariant of the redistribution process. 



3.5 Comparison with other work 

We calculated an atmosphere model for Tofi = 11 530 K, 
log q = 4.43, corre s pond ing to one of the envelope models of 



Landstreet et al, (1998). Then we calculated the radiative 
accelerations, but in the case of Nel we omitted all except 
the resonance lin es, to simulate the assumption adopted by 
This assumption is not really justified as 



Landstreet et al 



the omitted lines contribute significantly to the radiative 
acceleration. Inclusion of Ne I lines originating from excited 
levels, which absorb at the wavelengths where a star radiates 
much more than in the Lyman continuum, may increase the 
acceleration on Nel by more than 1.5 dex. However, this 
does not have any serious consequences on the conclusions of 
Landstreet et al. , as the radiative accelerations are so small 



(23) 



they will still remain below the gravity. The comparison is 
shown in Fig. ^. In this case we also considered Di — 1 
in Eq. ^ to follow their method. The agreement is quite 
good and small departures at the base of the atmosphere are 
mainly caused by slightly different temperature behaviour 
between our model atmosphere and their envelope model. 
The curious gap in our curve at 10~* < dm < 10~''gcm~^ 
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Landstreet et al. 



TjfpllSSO K, log g=4.43 



Budaj & Dworetsky ; 



10'" 1 
dm [cgs | 



10- 



Figure 8. Comparison of q u 



those of Landstreet et al. 



calc ulations in the atmosphere with 



(1998) in the envelope for the same 



model parameters. Here, for Nel only resonance lines are consid- 
ered. Both were calculated for standard abundances of Ne. 



appears because the acceleration is negative here and cannot 
be plotted on a logarithmic scale. 

3.6 NLTE effects 

The aforementioned negative acceleration is a very inter- 
esting NLTE effect which may lower still further the radia- 
tive accelerations on some ions in this part of the atmo- 
sphere. It is a consequence of the temperature inversion in 
the model atmosphere and the atomic structure of Nel. As 
shown in Fig. |l|, temperatures in the outer layers rise in a 
NLTE model for 10^** < dm < 10"^ gcm"^. At this depth, 
and in these stars, the atmosphere is still optically thick 
in the Lyman continuum. Thus, according to the diffusion 
approximation, the star radiates into itself towards cooler 
regions at these wavelengths. If an element happens to have 
its important transitions in the Lyman continuum it may 
be pushed into the star by the radiation and its accelera- 
tion is thus negative. Similar effects can be expected also in 
the cores of very strong lines not necessarily in the Lyman 
continuum. To explore NLTE effects on the results we also 
calculated radiative accelerations on Nel for three different 
assumptions in a representative NLTE model atmosphere 
with TcfT = 15 000 K, log p = 4 and zero microturbulence: 

(i) LTE case: neon populations in Eq. ^ are in LTE and 
the neon lines' source functions for flux calculations is equal 
to the Planck function 



(24) 



(ii) approximate NLTE: the same, but the line source 
function for resonance transitions is allowed to deviate from 
the Planck function following second order escape probabil- 
ity methods 



e+{l-e)K2{T\ 



(25) 



where e is the photon de struction probability and K2{t) is 



the kernel function (see fHubeny et al.| [1994 the SYNSPEC 



manual, for more detail). Note that this option was used in 
the calculations described in Section 2.2 above; 



Table 1. Nell energy levels considered. Column 1: level designa- 
tion, Column 2: ionization energy in c m~^. Column 3: statistical 



weight of the level. Energies are from Persson (1971) 
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Figure 9. Comparison of radiative accelerations on Nel obtained 
using different assumptions: LTE, approximate NLTE and full 
NLTE (see text). Calculated for standard abundance of Ne. 



(iii) fully consistent NLTE case: with NLTE model atmo- 
sphere, NLTE populations and NLTE source functions 



2hV^ 



gu,ini^, 
gi.iriu,: 



(26) 



In the fuU NLTE case the NLTE model atmosphere and Ne I 
level populations were calculated with tlusty195. As an in- 
put for TLUSTY195 we considered Hubeny's Hi atom model 
with 9 explicit lev els and continuum, t h e 31 level Ne I atom 
model taken from Dworetsky fc Budaj ( 2QQCl| ) , and a simple 
4 level atom of Ne 11 with continuum (Table ll|) . Generally, 
the data for Ne bound-bound and bound-free transitions 
were taken from the TOPbase ( CJunto et al.|p"993[ Hibbert 



Scott 1994) but Nel oscillator strengths were from Scaton 
(199S). The modion IDL interface written by Varosi et al 



( 1995 ) was particularly useful in constructing Ne I and Ne 11 
atom models. Populations of those very high Ne I levels not 
included in our Nel atom model were considered in LTE 
relative to the ground state of Nell, of which the popula- 
tion was calculated in NLTE, when calculating the energy 
flux and acceleration. All Nel level populations were set to 
LTE below dm > 0.3gcm~^ when calculating the radiative 
accelerations. 



The differences between the three assumptions are ex- 
hibited in Fig. ^ and become crucial for dm < lO^"* gcm"'^, 
where the radiative accelerations in LTE and NLTE may 
differ by more than 3 dex. In the case of LTE we observe 
negative acceleration in the region of the temperature in- 
version. This was explained above. The approximate NLTE 
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case is parallel to the LTE acceleration but about 1.5 dex 
smaller. This is caused by the difference in computed energy 
fluxes. While flux in the approximate NLTE case is real - 
lines are in absorption - and surprisingly close to the flux in 
the full NLTE case, the resonance Nel lines are in emission 
in LTE as soon as the atmosphere becomes optically thin in 
the Lyman continuum and acceleration soars immediately. 
The difference between approximate NLTE and full NLTE 
originates mainly in the differences in level populations. 



4 SUMMARY 



Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta 
53, 197 

Auer L.H., Mihalas D., 1973, ApJ 184, 151 
Babel J., 1995, A&A 301, 823 

Borsenberger J., Michaud G., Praderie F., 1979, A&A 76, 
287 

Budaj J., 1994, in: CP and Magnetic Stars, eds. Zverko J., 
Zizriovsky J., Astronomical Institute, Slovak Academy of 
Sciences, Tatranska Lomnica, Slovakia, p. 72 

Budaj J., Dworetsky M.M., Smalley B, 2002, Comm. Univ. 
London Obs. No. 82 



URL= http:/ /w ^ ww. ulo.ucl.ac.uk/ulo_comms 



We ha-i'e calculaLtid radlaiive acceleraliouH ou Nei-iv Iouh 
in the atmospheres of late B main sequence stars. We take 
into account the fine structure and calculations include the 
effects of line blending using the whole Kurucz line list, 
bound-free and free-free opacity of explicit elements (H, He, 
C) as well as some NLTE effects. We explored how the accel- 
eration changes with respect to effective temperature, sur- 
face gravity and homogeneous Ne abundance and found that 
it is much smaller than the gravity in and above the ob- 
served line forming region. Only at the base of our models 
(dm ~ 20gcm~^) could a Ne deficit of about 0.5-2.0 dex be 
supported by the radiation, depending on the effective tem- 
perature. This implies that, in the stable atmospheres of late 
B stars, Ne should be almost completely depleted from the 
photosphere. This is qualitatively i n agre ement with the ob- 
servations of Dworetsky & Budaj (2000) but the fact that 



/82/index.html 



we detected Ne in some HgMn stars suggests the presence 
of some weak transport or mixing mechanism contaminating 
observed layers by Ne from the reservoir underneath. Possi- 
bly, this mechanism is associated with the partial hydrogen 
ionization at the cool end and a very weak stellar wind at 
the hot end of the ffgMn T^s span. 

Finally, we demonstrated by a full NLTE calculation 
the importance of such effects on radiative accelerations and 
slightly improved a current treatment of the redistribution 
effect which should be used in future calculations. 
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